Filling-dependence of the zigzag Hubbard ladder for a quasi-one-dimensional 

superconductor Pr 2 Ba 4 Cu 7 Oi5_5 
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We investigate filling dependence of the zigzag Hubbard ladder, using density matrix renormal- 
ization group method. We illustrate the chemical-potential vs. electron-density and spin gap vs. 
electron density curves, which reflect characteristic properties of the electron state. On the basis of 
the obtained phase diagram, we discuss the connection to a novel quasi-one-dimensional supercon- 
ductor Pr2Ba4Cu70i5_« . 
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I. INTRODUCTION 

Low-dimensional strongly-correlated-electron systems 
have been one of the most active research fields in the 
condensed matter physics, since the discovery of high- 
ly superconductivity^ One of the current topics in the 
field is the frustration effect, which often causes a variety 
of interesting behavior such as exotic superconductivity. 
Recently, a novel quasi-one-dimensional(lD) compound 
Pr 2 Ba4Cu70i5_5(Pr247) was found to exhibit the super- 
conductivity below T c ~ 18K. This compound has the 
layered structure of single chains, frustrated zigzag lad- 
ders, and Cu02 planes^ An important point on Pr247 is 
that the metallic conductivity originates from the zigzag 
ladder part, while the Cu0 2 plane is insulating with the 
antiferromagnetic order below 285K i ' i 6 Moreover the 
oxygen atoms at the single chains, which are also in- 
sulating, are easily defected with the deoxidization and 
then the electrons are effectively doped into the zigzag 
ladders from the single chain sites. This suggests that 
Pr247 can be regarded as a natural "filling controlled 
system" , which hopefully provides an interesting physics 
cooperatively induced by the filling dependence and the 
frustration effect, fn fact, the transition temperature of 
Pr247 exhibits the systematical relation to the oxygen de- 
ficiency 8, implying that the electron-filling plays a cru- 
cial role for appearance of the superconductivity. 

A relevant model describing Pr247 is the zigzag Hub- 
bard ladder, whose Hamiltonian is given by 

H = y\tid irT Ci + i a + i 2 cJ<T c i+2cr + h.c] + U 2jnttnq,(l) 

in i 

where Ci a is the electron annihilation operator at i site 
with spin a and = c\ a Ci a ^ tifo) is the nearest(next- 
nearest) hopping term of the electron and U is the on site 
coulomb energy. We also introduce a — |ii/*2| for later 
convenience. According to Ref. 7 , the overlap of Cu-3d 
orbit and 0-2p orbit in the zigzag array of Cu-0 double 
chain structure suggests that the parameters correspond- 
ing to Pr247 are in t 2 < and a < 1. We also assume 
U = 8 1*2 1 or 8|ti|, which may be of similar order to the 
usual cuprates. The filling of the zigzag ladder part in 
Pr247 of 5 — corresponds to the nearly quarter filling. 5 



As the oxygen deficiency S increases, it rises continuously 
toward half-filled side and the superconductivity appears 
for 8 > 0.3. This indicates that it is primarily important 
to investigate the frustration and filling dependences of 
the electron state systematically, for through understand- 
ing of the Pr247 superconductivity. 

The Hamiltonian (JT]) is a typical example of the ID 
correlated electron system capturing various frustration 
effects. Indeed, interesting properties of low energy ex- 
citations are revealed at the half-filling or quarter fill- 
ing by extensive studies with bosonization and vari- 
ous numerical methods i&^ i 10 ' 11 : 12 ' 13 : 14 In particular, Fab- 
rizio sketched the qualitative phase diagram^ by invoking 
the weak coupling theory for the non-frustrating ladder 
systemi 5 ., and Daul etal presented an approximate phase 
diagram in the context of the ferromagnetism™, which 
also reveal that the zigzag Hubbard model contains quite 
rich physics. However, both the frustration and the in- 
commensurate Fermi wave number make a quantitative 
analysis of the low-energy excitations subtle; the precise 
filling-dependence of the spin gap is still unclear in the 
relevant parameter region to Pr247. 

In this paper, we precisely investigate the filling 
dependence of the zigzag Hubbard ladder with in- 
tensive calculations of density matrix renormalization 
group(DMRG) 16 . In the next section we particulary 
illustrate the chemical potential-electron density(/x-p) 
curve. In section III, we investigate the filling depen- 
dence of the spin gap in the region t% < 0, from which 
we read characteristic properties of the electron state . 
In sections IV and V, we summarize the DMRG results 
as a phase diagram and then discuss the relevance to the 
superconductivity in Pr247 respectively. Here, we note 
that the convergence of DMRG computation in the sin- 
gle band region is good with a relatively small number 
of retained bases m ~ 200, while in the two band region, 
we can obtain the reliable spin gap with up to m = 1000. 



II. fi-p CURVES 

Let us write the gourd state energy of L sites as 
El(N, S z ), where N denotes the number of electrons and 
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S z indicates the total-.!? 2 of the system. The electron 
density is written as p = N / L and then the chemical po- 
tential is defined as p = -[E L (N + 1, 1/2) - £ L (A,0)] 
for A^=even or /x = -[E L {N + 1,0) - El(N, 1/2)] for 
iV=odd. This definition of jit should be contrasted to the 
conventional notation of even number of electrons. This 
is because the p-p curve of the present definition of p can 
visualize some important natures of the low-energy exci- 
tation; If the charge excitation is gapless, the p-p curve 
acquires a smooth stairway-like curve. If two electrons 
conform a bound state, the N =even case is slightly sta- 
ble than the N =odd case, due to the binding energy of 
the electrons. This suggests that if the system has the 
bound state, the p-p curve exhibits overhung behavior. 
We should note that the similar behavior can be actu- 
ally seen in the magnetization curve of the zigzag spin 
system.— 

As can be expected, the basic property of the p-p curve 
is attributed to the one-particle dispersion curve that is 
easily obtained as 



e{k) = 2<i cos(fc) + 2i 2 cos(2fc). (2) 



The shape of this dispersion curve has double-well form 
for a < 4 (recall the sign of t% is negative), where it has 
the van Hove singularity. In the following, the corre- 
sponding electron density is denoted as p c . As a in- 
creases, p c comes down to the lower-electron density re- 
gion. For a < 2, p c reaches below the half filling^ < 1), 
and, at a = 1/2, it is located near the quarter fill- 
ing. As a — > 0, the system finally becomes the two 
decoupled Hubbard chains. Thus the system in the non- 
interacting limit (J7 = 0) is essentially a single band sys- 
tem for p < p c , while it is a two band system having four 
Fermi points for p > p c . We note that a finite U may 
induce the spin gap in the two band region, according to 
the weak coupling theory. 

In Fig. [TJ we show the p-p curves of U — 8 for a = 1 
and a = 0.5 with t% = — 1. In the figure, we can find 
a characteristic feature: The curvature of the p-p curve 
changes at p c ~ 0.73 for a = 1 and p c ~ 0.5 for a = 0.5, 
which is basically reflecting the shape of the dispersion 
curve We however note that these values of p c are 
slightly lifted from the non-interacting case due to the 
effect of U. A more interesting point is that the p-p 
curves actually exhibits the oscillating behavior above 
p c . For a = 1, we can see the overhung behavior in 
p c < p < 1, implying that the spin gap may exist between 
p c and the half filling. On the other hand, for a = 0.5, 
the oscillation is enhanced above p c , but the p-p curve 
becomes smooth again above p ~ 0.83 as is indicated by 
an arrow in FigQ] (b) . This suggests that the spin gap for 
a = 0.5 may appear only near p c . Here, we should note 
that the even-odd effect below p c is properly removed by 
the size extrapolation. 




FIG. 1: p-p curves for L = 60 systems of U = 8: (a) ti = 1 
and t 2 = -l(a = 1) and (b) h = 0.5 and t 2 = -l(a = 0.5). 
The broken lines are guides for p c . 

III. SPIN GAP 

Since we have seen the outline of the filling dependence 
of the system in the p-p curve, we further analyze the 
spin gap for the precise characterization of the low-energy 
excitation. We define the spin gap A s as 

A S (N;L) = E L (N,1) ~ E L (N,0), (3) 

for N =even. In order to extract the bulk behavior of 
the spin gap, the finite size extrapolation is usually re- 
quired. We should however recall that p takes only some 
fractional values restricted by combination of L and N, 
as far as we treat a finite size system. In order to avoid 
this mismatching problem of the electron filling, we cal- 
culate A s of all the electron numbers for various L, with 
which we interpolate the spin gap at an irrational value 
of p. We think that the bulk spin gap can be properly 
analyzed except for the vicinities of some singular points. 

In Fig. [2] (a) , we show the DMRG calculated spin gap 
for L = 48,60,72 with t x = -t 2 = l(a = 1). We can 
see that the property of the spin gap clearly changes at 
p c ~ 0.7, which is consistent with the p-p curve. Be- 
low pc, the system is essentially described by the single 
band and A s shows large size dependence, suggesting 
that the spin excitation is gapless. We thus examine 
the 1/L-size dependence of the interpolated spin gap: 
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FIG. 2: Size dependence of spin gaps for U = 8: (a) ti = 1 
and t2 = —1 and (b) ti — 0.5 and £2 = — 1- 



A s (p;L) — A s (p) + const/L, and then verify that A s 
below p c becomes zero in L — > oo. On the other hand, 
we can expect the spin gap in p c < p < 1, where the 
overhung behavior of the p-p curves is observed. In fact 
we can clearly see that A s is almost independent of the 
system size, implying that the spin gap remains a finite 
value in the bulk limit. Moreover, we can find that the 
spin gap is enhanced in p < 1 rather than the half-filling, 
which may be a peculiar behavior in the frustrating sys- 
tem in contrast to the non-frustrating Hubbard ladder. 
Here, we note that this spin gap phase is adiabatically 
connected to the dimer spin gap of the corresponding 
zigzag spin system at the half filling. The amplitude of 
the spin gap at the half filling is consistent with the zigzag 
spin chain 18 . 

Fig. [U (b) shows A s for L = 48, 60, 72 with t x = 0.5 
and ti — — 1(q = 0.5). For p < p c {— 0.5), the system is 
essentially a single band model and thus the size depen- 
dence of A s is basically the same as that for a = 1; the 
1/L extrapolation of A s leads that the spin excitation is 
gapless in the thermodynamic limit. On the other hand, 
A s above p c shows subtle behaviors. Here, we should re- 
call that the p-p curve for a = 0.5 has a weak anomaly 
at p ~ 0.83, above which p-p becomes smooth. We first 
analyze the region of p c < p < 0.83, for which the even- 
odd oscillation appears in the p-p curve. In this region, 
the spin gap shows oscillating behavior with respect to 
p and a. The amplitude of the oscillation increases, as 



p increases from p c to 0.83. Thus the precise extrapo- 
lation of A s is still difficult in this region. However, we 
note that the size dependence of A s in the vicinity of 
Pc{Pc < P < 0.7) is rather weak and a small spin gap 
may survive in the bulk limit. 

Next let us turn to 0.83 < p < 1, where we can see 
the similar oscillating behavior of A s . However, this os- 
cillation shows rather systematical behavior depending 
on N mod 4 = or ^ 0. In this region, the system is 
basically described by the two chains, where the electron 
filling is sufficiently away from p c . Since the system has 
four Fermi points, the N mod 4 = state is stabilized by 
conforming "closed shell" with respect to the four Fermi 
points, while the N mod 4^0 state has a "unoccupied 
orbits" , which may generate an anomalous spin excita- 
tion. We thus read the bulk behavior of the spin gap from 
N mod 4 = sectors, which shows large size dependence 
suggesting that the spin excitation seems to become gap- 
less or very small. However, the precise estimation of the 
spin gap value is also difficult within L — 72. 

Here, we make a comment on the boundary effect, since 
a boundary excitation may appear in the gapped frustrat- 
ing system with the open boundary. In order to check it, 
we have also calculated magnetization curve at some typi- 
cal fillings. Then we can verify that the S z = 1 excitation 
smoothly connects to the bulk part of the magnetization 
curve in the spin gap phase of a = 1 . On the other hand, 
the magnetization curve of a — 0.5 suggests that the A s 
sometimes capture the boundary effect in the oscillating 
spin gap region, where the S z — 1 excitation sometimes 
generates an anomalous step away from the bulk part of 
the magnetization curve. 



IV. PHASE DIAGRAM 

In Fig. [3l we summarize the phase diagram of the 
zigzag Hubbard model for U — 8, which is determined 
by the p-p curves and the spin gap for L = 72. The 
open circles indicate the boundary between the single- 
band gapless phase and spin gap/oscillating region, which 
is slightly lifted from p c for the free electron case due to 
the electron correlation effect. The " spin gap" indicates 
the spin gap phase, which is connected to the dimer spin 
gap phase at the half filling^^. As a decreases from 
the dimer spin gap phase, the spin gap rapidly decreases 
and almost vanishes across the dotted line. We also note 
that spin gap at the half-filling for a < 0.7 seems to be 
almost gapless. However, it is still difficult to distinguish 
whether the spin excitation in this regime is the truly 
gapless phase or very small spin gap phase within the 
present accuracy. 

In the region "oscillating" , the behaviors of the spin 
gap becomes subtle, as in FIGlDJb). The oscillating be- 
havior of A s makes precise analysis of the spin gap in 
the bulk limit difficult. In the vicinity of p c , however, we 
have seen that the size dependence of A s is not so large, 
suggesting the existence of the spin gap, which is relevant 
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FIG. 3: Spin-excitation phase diagram of the zigzag Hubbard 
ladder of U = 8. In the left-half, < -ti < 1 with ti = 1, 
and in the right half, < ti < 1 with ti = — 1. Open circles 
indicate the boundary between the single-band gapless phase 
and the spin gap phase. The dotted line is a guide for the 
vanishing line of the spin gap. Triangles mean the boundary 
between the (almost)gapless region and the oscillating spin 
gap region. The broken line means p c of the free electron 
case. The error bar corresponds to the interval of descritized 
density of electron for a L = 72 system. 



to Pr247(see Sec. V). As p increases, the irregular oscil- 
lation to p becomes significant and the size extrapolation 
of A s breaks down. The upper bound of the oscillat- 
ing region is indicated by the open triangles, which are 
corresponding to the weak kink in the p-p curve. In or- 
der to illustrate another aspect of the oscillating region, 
we further investigate A s as a function of a for a fixed 
p. Figure |4] shows the a-dependence of A s for L = 72 
and N = 60{p = 0.833). We can then find that A s ex- 
hibits the rapid oscillation with respect to a for a < 0.5, 
which suggests that the incommensurate nature due to 
the frustration and the boundary effect cooperatively in- 
duce such subtle behavior of A s . The rapid increase of 
A s for a > 0.8 corresponds to the spin gap phase in Fig. 
[3l In the intermediate region(0.5 < a < 0.7) the size 
dependence of A s suggests gapless or almost gapless spin 
excitation, as mentioned before. 

Here we would like to make a comment on the res- 
olution of the phase diagram within L = 72; For ex- 
ample we can not distinguish whether the boundary be- 
tween the gapless region and the oscillating region/spin 
gap phase(open triangles/dotted line in Fig[3]) is a bulk 
phase transition or not, chiefly because of the limited sys- 
tem size. Also such a narrow phase as C2S2 mentioned 
by the weak coupling theory 8 is beyond the resolution of 
the discretized density of the electron. In addition, we 
note that the spin gap for a < 0.3 is not evaluated prop- 
erly within the DMRG calculation, since the energy scale 
of the spin gap itself becomes too small in the decoupled 
chains limit. 
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FIG. 4: a-dependence of the spin gap A 3 of the zigzag Hub- 
bard model for L = 72 and TV = 60(p = 0.833). The 
on-site Coulomb interaction is fixed at U = 8. In a > 
0.7(corresponding to the dotted line in FIG[3]), the system 
is in the spin gap phase. In a < 0.5, A s shows the rapid 
oscillation with respect to a. 



V. DISCUSSIONS 

On the basis of the phase diagram, let us discuss the 
filling dependence of Pr247. As was discussed in Ref. 7 , 
the chemical potential shifts from the quarter filling to- 
ward the electron doping side by oxygen reduction. If 
assuming that T c is proportional to the amplitude of the 
spin gap in the zigzag ladder, we can see that the effective 
model parameter of Pr247 is most likely in a ~ 0.5 ~ 0.7, 
for which the S dependence of T c basically agrees with 
the filling dependence of the spin gap 3 . Thus, Pr247 
may be located near the most competing region in the 
phase diagram where the spin gap shows the subtle be- 
havior, which is consistent with the relatively low T c 
of Pr247. In addition, it is suggested that Pr247 is 
sensitive to modification of parameters due to external 
fields such as high pressure effect^ As confirmed in the 
NQR/NMR experiments^, a good one-dimensionality can 
be expected for Pr247 due to the superstructure of the 
metallic zigzag ladders, the insulating single chains and 
Cu02 planes. Such superstructure might be related to 
the relatively large a, compared with a similar ladder 
compound YBa2Cu40s- Very recently, Nakano et aZhave 
made FLEX calculation on the basis of the band struc- 
ture of YP^CujOs, which nevertheless suggests the s- 
wave superconductivity in the small a region^. Thus 
the consistency with the actual band structure of Pr247, 
which is not available now experimentally and theoreti- 
cally, is highly interesting. 

In this paper, we have clarified the various interesting 
behaviors induced by the frustration effect in the phase 
diagram. In particular, the enhancement of the spin gap 
above p c and the appearance of the oscillating region 
may be essential in the connection with the Pr247 ex- 
periments. However, we did not mention the correlation 
functions and the pairing symmetry here, which has been 
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still unknown experimentally. In order to determine the 
paring symmetry, it is needed to precisely investigate the 
correlation functions for a sufficient long chain. Also the 
Independence of the zigzag Hubbard model, e.g. the con- 
nections to the t- J model and the weak coupling theory, is 
a theoretically important problem. We hope the present 
work to be a portal of further researches on the Pr247 su- 
perconductivity from both theoretical and experimental 
view points. 
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